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ABSTRACT 

We present J, H and K-band spectroscopy of Cygnus A, spanning 1.0-2.4/zm in the 
rest- frame and hence several rovibrational H 2 , H recombination and [Fell] emission 
lines. The lines are spatially extended by up to 6kpc from the nucleus, but their 
distinct kinematics indicate that the three groups (H, H2 and [Fell]) are not wholly 
produced in the same gas. The broadest line, [FeII]A1.644, exhibits a non-gaussian 
profile with a broad base (FWHM ~ 1040 kms -1 ), perhaps due to interaction with 
the radio source. Extinctions to the line-emitting regions substantially exceed earlier 
measurements based on optical H recombination lines. 

Hard X-rays from the quasar nucleus are likely to dominate the excitation of the 
H 2 emission. The results of Maloney, Hollenbach & Tielens (1996) are thus used to infer 
the total mass of gas in H2 v=l-0 S(l)-emitting clouds as a function of radius, for gas 
densities of 10 3 and 10 5 cm -3 , and stopping column densities TVh = 10 22 — 10 24 cm~ 2 . 
Assuming azimuthal symmetry, at least 2.3 x 10 s M Q of such material is present within 
5 kpc of the nucleus, if the line-emitting clouds see an unobscured quasar spectrum. 
Alternatively, if the bulk of the X-ray absorption to the nucleus inferred by Ueno 
et al. (1994) actually arises in a circumnuclear torus, the implied gas mass rises to 
~ 10 10 M Q . The latter plausibly accounts for 10 9 yr of mass deposition from the cluster 
cooling flow, for which M ~ 10 M Q yr^ 1 within this radius. 
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1 INTRODUCTION 

X-ray observations show that the gas within the inner few 
100 kpc of many clusters of galaxies is cooling from tem- 
peratures of 10 7 — 10 8 K at rates of up to several hun- 
dred solar masses per year. If persistent for 10 s yr or more, 
the inferred subsonic inflow of cooling material- the cool- 
ing flow (CF) (see Fabian 1994 for a review)- would deposit 
in excess of 10 10 Mq of cool gas around the central cluster 
galaxy (CCG). Whilst the presence of such large amounts 
of cold material has been deduced from soft X-ray absorp- 
tion (White et al. 1991, Allen et al. 1993, Allen et al. 1995, 
Allen & Fabian 1997), its detection in emission has proved 
problematic and the ultimate fate of the gas remains uncer- 
tain. The bulk of it is clearly not turned into visible stars, 
but the role of dust in obscuring any star formation may 
be important (Allen 1995; Crawford et al. 1999), and recent 
sub-millimetre work in fact suggests star formation rates of 
up to lOOMoyr -1 (Edge et al. 1999). 

There is also a strong link between the cooling time 
in the central 30 kpc of a cluster and the presence of optical 



line-emission and/or a CCG radio source: Peres et al. (1998) 
show that for central cooling times below 3 x 10 9 yr, 90 per 
cent of the CCGs exhibit line emission, and 95 per cent 
contain a radio source; for comparison, in CFs with longer 
central cooling times, the fractions are 5 and 20 per cent, re- 
spectively. Despite the uncertainty over the mechanism be- 
hind the optical line emission, the surrounding ICM clearly 
plays an important role. There have recently been attempts 
to study the emission from cooler phases of the gas, such as 
molecular hydrogen which emits through thermal excitation 
at ~ 2000 K. Jaffe and Bremer (1997) detected H 2 v=l- 
S(l) emission in the inner few kpc of three CF CCGs, 
but not in any radio galaxies outside strong CFs at simi- 
lar radio powers and redshifts. They found that the emis- 
sion is too luminous to be due to material simply passing 
through ~ 2000 K whilst cooling from higher temperatures, 
and implicated suprathermal electrons deep in the molecu- 
lar core of the cloud as the reheating mechanism. Falcke et 
al. (1998) also observed strong H2 emission associated with 
radio galaxies in CFs, and speculated that it is due to radio- 
jet interactions with the ICM or to the injection of material 
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from infalling spiral galaxies. Krabbe et al. (2000) presented 
H and K-band spectral imaging of NGC 1275, the cD galaxy 
at the centre of the Perseus cluster, a 200Moyr -1 CF. 
They cited the central concentration of the H2 emission, its 
relatively low and uniform temperature (1500-3000 K) and 
AGN-related excitation mechanism as evidence against the 
CF hypothesis. But such properties are in fact entirely con- 
sistent with the expectation that much of the gas is too cold 
(< 100 K) to emit significant H2 emission, and is only ob- 
served when excited in some way, in this case by the AGN 
itself. 

We have recently begun a near-infrared study of the 
H2 emission in a sample of 28 CF CCGs, chosen to be the 
most optically line luminous in the complete spectral study 
of 217 ROSAT-selected clusters of Crawford et al. (1999). In 
this paper we present J, H and K band spectra of Cygnus 
A, which was observed early in the programme. At redshift 
z = 0.056, it is by some considerable margin the most pow- 
erful radio source in the local universe and, like the radio 
galaxy 3C295, similar in power to the z > 1 FRII radio 
sources. It resides in a moderately rich, cooling flow cluster 
with a mass deposition rate of ~ 250 Mq yr _1 (Reynolds & 
Fabian 1996); X-ray observations also reveal an absorbed 
power-law component from the quasar nucleus (Ueno et 
al. 1994), which further manifests itself through scattered 
broad emission lines (Antonucci, Hurt & Kinney 1994; Ogle 
et al. 1997). 

This paper is structured as follows: the observations and 
data reduction are described in section 2, followed in sec- 
tion 3 by a discussion of the extended emission in the lines 
of molecular hydrogen, the hydrogen recombination lines of 
Paa, Pa7, Br7 and Br<5, and [Fell] AA1. 258, 1.644; in section 
4, we make deductions about the gas content and molec- 
ular excitation mechanisms within the central few kpc of 
Cygnus A, in the context of the current understanding of the 
CF and the obscured quasar. The cosmological parameters 
Hq — 50kms _1 Mpc -1 and qo — 0.5 are adopted through- 
out, yielding a spatial scale of 1.5 kpc arcsec - at the redshift 
of Cygnus A. 



2 OBSERVATIONS AND DATA REDUCTION 

The observations were taken with the CGS4 spectrograph 
on the United Kingdom Infrared Telescope (UKIRT) in 
September and October 1999, as shown in the observation 
log in Table 1. The 256 x 256 InSb array, 40 1/mm grat- 
ing and 300 mm focal length camera were in place, yielding 
a spatial scale of 0.61 arcsec per pixel and spectral resolu- 
tions of 610, 880 and 570 kms" 1 FWHM (in the J, H and K 
bands, respectively) with a 2-pixel wide slit aligned north- 
south. The NDSTARE mode was used along with the con- 
ventional object-sky-sky-object nodding pattern, thus ob- 
viating the need for separate bias and dark current frames 
and permitting the computation of an external error for each 
pixel; it also reduces the required accuracy of the flat-field 
frame. Atmospheric absorption features were removed by 
ratioing with some of the main sequence F stars tabulated 
at the telescope which were calibrated against photometric 
standards (except for the J-band exposure acquired under 
non-photometric conditions, and for which an approximate 
flux calibration was obtained by matching the continuum 



level with that in the H-band). On- and off-line data re- 
duction was performed using version VI. 3-0 of the Portable 
CGS4 Data Reduction package available through Starlink. 
Row-by-row spectra were extracted from the fully-reduced 
spectral images and converted to ASCII format for use with 
the emission-line fitting package QDP/PLT (Tennant 1991). 
Emission lines were fitted with gaussian components atop a 
linear continuum. 



3 RESULTS 

Spectra from the nuclear row of Cygnus A are shown in 
Fig. |l|. Many of the lines were previously identified by 
Ward et al. (1991) and Thornton et al. (1999) in ear- 
lier UKIRT spectra. The blue shoulder on Paa is plau- 
sibly due to HeIIAl.8639 (using case B theory and the 
HeIIA4686 flux in Osterbrock & Miller 1975, its pre- 
dicted flux is 3 x 10~ 16 erg cm~ 2 s" 1 ) and also perhaps 
to HeIAAl.8691, 1.8702 (although we cannot predict their 
fluxes); Thornton et al. (1999), however, ascribed it en- 
tirely to very high excitation H2 O-series lines. Using the 
list of coronal lines in Ferguson et al. (1997), we identify 
several new such lines, viz. [Ca VIII]A2.32, [Si XJA1.43 and 
[S IXJA1.252, in addition to those previously known between 

2.0 and 2.1/zm. As an alternative to [S IXJA1.252, the blue 
component of the doublet near 1.325/^m could instead be 
HeIAl.2531 and/or H 2 v=2-0 Q(4), but if it were the latter, 
the absence of the Q(l), Q(3) and Q(5) lines of this series 
(which should be much stronger under thermal excitation) 
would be hard to explain. In the present data, many of the 
identified emission lines extend significantly beyond the in- 
strumental point spread function (PSF), as described below. 

3.1 Spatially extended line fluxes 

Fig. ^ shows the emission line fluxes as functions of posi- 
tion along the slit. The nuclear light of the H-band expo- 
sure falls on two rows of the chip, so the [FeII]A1.644 points 
in Figs. |^ and ^ have been shifted by +0.5 pixel. All the 
lines are spatially extended, especially to the north, and by 
up to 6.3 kpc in the case of v=l-0 S(l). There are notable 
variations in the line ratios, e.g. in S(l)/S(3), and in ratios 
involving [Fell] AA1.258, 1.644, Br 7 and v=l-0 S(l), which 
we discuss in section 4. We caution, however, that there are 
systematic uncertainties in the S(3) flux in rows 96-98 owing 
to the blending of this line with [Si VI]A1.962. 

3.2 Kinematics 

Fig. ^ shows the variations along the slit in the radial ve- 
locity and FWHM of H 2 v=l-0 S(l), Paa and [Fell] A1.644. 
That the kinematics of the three lines differ from one an- 
other shows that they are not entirely produced in the same 
gas. The radial velocity curves of H2 v=l-0 S(l) and Paa, 
showing a rise to the north to ~ 400kms _1 and an initial 
decrease to the south before increasing again, are at least 
qualitatively consistent with the behaviour of [OIII]A5007 
and [NII]A6584 found by Stockton et al. (1994) at the same 
PA. Using a narrow slit offset from the nucleus by various 
amounts, they found that these lines split into two compo- 
nents separated by 300kms _1 1 arcsec west of the nucleus 
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Table 1. Log of observations 



UT Date A — A Resolution f SeeingJ Integration Slit 

(1999) (fan) (kms -1 at FWHM) (arcsec) time (min) PA (deg) 

Sep 2 1.875-2.490 (K) 570 1.0 64 

Sep 3 1.345-1.961 (H) 880 1.2 48 

Oct 7 1.128-1.435 (J) 610 1.5 * 52 

f At centre of the wavelength range. 

| FWHM of standard star spatial profile along the slit. 
* The night was non-photometric. 



and that the gradient of the velocity curve also changes sign 
there, due to an emission component located north-west of 
the nucleus. Also from optical lines, Tadhunter et al. (1994) 
found relatively large line widths (200-600 kms" 1 FWHM) 
along the radio axis (PA 105 degrees), but only small radial 
velocity shifts (< 200 kms" 1 ), although Tadhunter (1991) 
had earlier found a high velocity cloud, moving at 1500- 
1800 kms" 1 from the galaxy rest frame and associated with 
the north-west emission line component. 

The S(l) line has FWHM~ 420 kms" 1 on nucleus (in 
agreement with the measurement by Thornton et al. 1999), 
and is significantly broader to the south than to the north. 
In rows 96-98 where S(3) is strongly blended with [Si 
VI]A1.962, the velocity width of the former was set to that 
of the S(l) line in a two gaussian fit; the fitted width of [Si 
VI] was found to be ~ 930 kms" 1 FWHM in rows 96-98 
(consistent with Ward et al.'s value of 1200 ± 300 kms" 1 ), 
~ 550 kms" 1 in row 95 and unresolved further to the north. 
Thornton et al.'s low value of 582 kms" 1 for the FWHM 
of [SiVI] on nucleus and high value of 552 kms" 1 for the 
S(3) line (much higher than the other v=l-0 S lines) demon- 
strates that leaving both line widths free in fitting such a 
blend can lead to errors (their value of 1.2 for the v=l- 
S(3)/S(l) line ratio is thus too high). We also note 
that on nucleus (row 96), the para- hydrogen lines v=l- 
S(0), S(2) and S(4) are all significantly narrower than 
the ortho-hydrogen lines, v=l-0 S(l), S(3) and S(5): S(l), 
S(3) and S(5) have a mean width of 410 kms" 1 FWHM, 
whereas S(0) and S(2) are unresolved and the S(4) width 
~ 200 kms" 1 FWHM; this implies that in this row the ortho 
and para lines may in part be produced in distinct regions 
with different (and at least in part, non-thermal) ortho/para 
population ratios. 

The fitted width of Paa is ~ 630 kms" 1 FWHM on 
nucleus, and generally higher to the south than to the north 
(its blue shoulder was modelled with an additional compo- 
nent). The [FeII]A1.644 line is as broad as Paa on nucleus 
but shows an off-nucleus increase, which may be due to an 
additional narrow component on nucleus. This is shown in 
Fig. § for row 95, where the broad (FWHM=1040kms" 1 ) 
base contains 70 per cent of the flux, with the narrow (un- 
resolved) core providing the rest. 

We speculate that the broader, spatially-extended 
[FeII]A1.644 component results from shocks (e.g. from the 
radio source), which plausibly produce the larger line- width 
and raise the gas phase Fe abundance by evapouring it 
from dust grains. In this regard, we note that Forbes & 
Ward (1993) attributed the [Fell] emission in starburst and 
Seyfert galaxies to shock-excitation in supernova remnants 
and radio jets, respectively. Like the latter authors, Simp- 



son et al. (1996) also found that the [Fell] emission in ac- 
tive galaxies is more tightly correlated with radio power 
than is the H recombination emission, but concluded that 
photoionization is the dominant [Fell] production mecha- 
nism, although ~ 20 per cent of it may be shock-excited. 
Krabbe et al. (2000) also observed a markedly non-gaussian 
[FeII]A1.644 profile in NGC 1275, leading them to suggest 
that it is produced in two distinct regions. We also note that 
Whittle (1992) found that the FWHM of the [OHI]A5007 
line in Seyfert galaxies without linear radio sources is largely 
determined by gravity (correlating well with the amplitude 
of the galaxy rotation curve); Seyferts with luminous linear 
radio sources, however, are observed to have significantly 
broader lines, due to acceleration by the radio source. 

In the J-band, the FWHM of Pa 7 
is 430^43Q kms" 1 (with la errors), although with the er- 
rors it is consistent with the value of 630 kms" 1 FWHM 
observed for Paa. The fitting of the [FeII]A1.258 line is com- 
plicated by the line on its blue wing which, as mentioned 
earlier, is either HeIAl.2531 or [SIXJA1.252: on nucleus, the 
fitted FWHM is 430^^ kms" 1 if the widths of the two 
lines are left free, or 570±H° kms 1 if they are constrained 
to be the same. These ranges are consistent with the FWHM 
of [FeII]A1.644 on nucleus, but it is not possible from these 
data to infer whether or not [FeII]A1.258 also has a broad 
base (since the two lines originate from the same upper level, 
they should have identical profiles in the absence of any dif- 
ferential extinction - recall from Table 1, however, that the 
instrumental resolution is slightly poorer in the H-band). 
This is illustrated by the line profile comparison in Fig. Be 
the lines have been assumed to have the same velocity cen- 
troid, although the fitted position of [FeII]A1.258 appears 
to be blueshifted by ~ 200kms" 1 from [FeII]A1.644 (but 
the HeIAl.2531 and/or [SIXJA1.252 would tend to reduce 
the fitted centroid of [FeII]A1.258, and the relative wave- 
length calibration between 1.25 and 1.64/im is only accu- 
rate to ~ 100 kms" 1 ). It should also be noted that the slit 
could have been in a slightly different position during the J- 
band observation, which was taken under non-photometric 
conditions with poorer seeing. 



3.3 Reddening 

On nucleus, we measure Br7/Paa = 0.107 ± 0.007 which 
translates into an extinction of Ay — 9.6 mag, assuming an 
intrinsic line ratio of 0.083 (case B) and the near-infrared in- 
terstellar extinction curve of Landini et al. (1984). For com- 
parison, Ward et al. (1991) deduced Av — 1.2 ± 0.15 mag 
using Paa/H/3, and Tadhunter et al. (1994) found Av — 3 
from Ha/H/3. There are two possible reasons for the dis- 
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[Fell] AM. 258. 1.644, v=l-0 S(l) and Bry line fluxes in Cygnus A 
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Figure 1. J, H and K-band spectra of the nucleus of Cygnus Figure 2. Emission line fluxes as a function of row number; 
A, plotted with la error bars obtained from the data re- the peak in the underlying continuum in row 96 of the In- 
duction; Pa/3 lies in the atmospheric absorption band at band exposure is taken to mark the position of the nucleus 
1.354/im. (the H-band points, i.e. [Fell] Al. 644, have thus been shifted 

by 0.5 pixel, as explained in the text; e.g. the line at row 96 
on the chip is shown here at 96.5). The unlabelled curves in 
the upper plot are PSFs in the intervals 1.9-1.95/im (dotted 
line) and 2.15-2.25/im (dashed line). The spatial scale is 0.6 
arcsec per pixel. 
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(a) Radial velocity relative to nuclear v=l— S(l) 



[Fell] A1.644 




95 

Row number 



(b) FWHM 



[Fell] M.844 



v=l-0 S(l)/\ 

— / 




90 95 100 

N Row number S 

Figure 3. (a) The radial velocities of the Hb v=l-0 S(l), Paa 
and [FeII]A1.644 as functions of row number, with respect 
to the v=l-0 S(l) emission in the nucleus (row 96). (b) The 
FWHM of these lines, after correction for the instrumental 
resolution (points at FWHM=0 kms -1 and without error 
bars represent unresolved lines). The spatial scale is 0.6 arc- 
sec per pixel. 



crepancies: firstly, Br7 may be contaminated by some Hel 
lines on its blue wing (although their expected strengths are 
not known); alternatively, the obscuring material may not 
be in the form of a uniform foreground screen (e.g. it may 
have variable thickness, or be mixed with the gas), in which 
case decrements between different optical and near-infrared 
lines would not yield the same Av- For completeness, we 
note that Br7/Paa = 0.112 ± 0.013 in row 94, equivalent 
to A v ~ 11.5 mag; in row 98, Br7/Paa = 0.0685 ± 0.033, 
consistent with the case B value and therefore with no ex- 
tinction. 

Drawing upon the J-band data, we measure Pa7/Paa = 
0.083 ± 0.019 on nucleus, which translates into Av = 
5.7 mag; in rows 95 and 97 we infer Av = 8.2 and 
5.4 mag, respectively. We may also compute the extinc- 
tion to the [FeII]-emitting region. Intrinsically, we expect 
a fixed [FeII]A1.258/A1.644 ratio of 1.35 (Nussbaumer & 
Storey 1988), but find on nucleus that [FeII]A1.258/A1.644 ~ 




1.7 1.72 1.74 1.76 

Observed wavelength (jam) 

Figure 4. The distinctly non-gaussian profile of the 
[Fell] A1.6441ine in row 95 of the H-band exposure; in 
a two component fit, the broader component (FWHM~ 
1040 kms~ ) contains 70 per cent of the flux, the rest being 
in the unresolved core. 




fab . :n □ 



AV (km s ') 

Figure 5. Profiles of the [FeII]A1.258 (solid) and [FeII]A1.644 
(dotted) lines in velocity space, after subtraction of the fit- 
ted linear continuum and scaled to the peak flux of the fitted 
line. The velocity centroid of each line is taken from a gaus- 
sian fit. It is difficult to infer whether the broad base to 
[FeII]A1.644 is also present on [FeII]A1.258, due to the pres- 
ence of either [S IXJA1.252 and/or HeIAl.2531 on the blue 
wing of [FeII]A1.258. 



0.3, implying Av = 12 mag; the extinction is at least as large 
in row 97, but somewhat less in row 95 (Av — 9 mag). As 
noted at the end of section 3.2, however, the presence of the 
broad component to [FeII]A1.644 which may be lacking from 
[Fell] Al. 258, suggests that this component is subject to more 
extinction than the values just calculated, and the narrow 
component which dominates [FeII]A1.258 to less (i.e. a sin- 
gle foreground extinction screen is not a satisfactory model) . 
We also reiterate the caveat concerning the possibility that 
the slit may have been at a slightly different position during 
the J-band observation. 

Concerning the extinction to the dust continuum- 
emitting regions, Imanishi & Ueno (2000) measured the opti- 
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Table 2. Extinctions to line-emitting regions 
Line ratio reference! 



Br7/Pao 
Pa7 /Paa 

Paa/H/3 Ward et al. (1991) 

Hq/H/3 Tadhunter et al. (1994) 

FeIIA1.258/A1.644 



LTE temperature diagnostic diagram for C)^gnus A 



A v 

on nucleus (mag) $ 

9.6 
5.7 
1.2 

3 
12 



f This work unless stated. 

X Calculated using a foreground screen extinction model. 



cal depth of the 9.7/im silicate absorption feature and found 
i"9.7 ~ 1 (equivalent to Av ~ 9 — 19 mag if the dust prop- 
erties resemble those of Galactic dust). That this figure is 
much less than the Av ~ 150 mag to the 3/im-emitting re- 
gion implies, they argue, that there is a dust temperature 
gradient, and their models suggest that this requires a dusty 
torus with an inner region of < 10 pc. It is of interest that 
their Av ~ 9 — 19 mag to the 10/im continuum emission site 
does not significantly exceed our inferred extinctions. For 
convenience, we summarise the latter in Table 2. 



4 INTERPRETATION 
4.1 LTE Comparison 

From the findings of section 3 we wish to infer the prop- 
erties of the line-emitting H2 and the associated excita- 
tion mechanism. If the v=l-0 S(l) emission is thermally ex- 
cited at a vibrational temperature of 2000 K, we may follow 
Scoville et al. (1982) and deduce the mass of warm H2 at 
this temperature; using the total v=l-0 S(l) luminosity of 
(5.7 ± 0.3) x 10 40 erg s" 1 (integrated along the slit), we find 
M(warm H 2 ) = (2.4 ± 0.1) x 1O 5 M , much less than the 
total mass which the cooling flow would have deposited, but 
the bulk of that material would in any case be much too 
cold (< 100 K) to produce thermally-excited S(l) emission. 
Following Jaffe & Bremer (1997), we may also compute the 
mass cooling rate, M, by multiplying the v=l-0 S(l) lu- 
minosity in photon s _1 by the mass of the H2 molecule and 
dividing by the fraction of the total cooling rate at 2000 K in 
this line, usually estimated at between 2 and 10 per cent: we 
find M > 3 x 10 4 M yr _1 , over 100 times greater than the 
X-ray mass deposition rate. We are thus driven to the same 
conclusions as Jaffe & Bremer, namely that the observed 
line emission is not simply due to material passing through 
~ 2000 K whilst cooling from 10 7 K and that an additional 
heat source is therefore required. 

Clues to the nature of the latter may be gleaned from 
Fig. ^, which shows for each row flux ratios involving the 
v=l-0 S(l), S(3) and S(5) lines; the upper levels of these 
transitions lie 6956, 8365 and 10341 K above the ground 
state, respectively, giving sensitivity to a broad range of 
temperatures. It is apparent from this figure that much of 
the emission is not thermally excited in high density gas 
at a single temperature, especially in the nucleus and to 
the south of it. The points deviate from the LTE predic- 
tion in the sense expected for emission from gas at a range 




0.5 



1 

ratio of v=f-0 S(3)/S(l) 



1.5 



Figure 6. Comparison of the ratios of several H2 v=l-0 S- 
series lines with the predictions for thermal excitation at 
various temperatures; the star denotes the nucleus (row 96) , 
circles rows to the south (rows 97 and 98) , and squares rows 
to the north (rows 92-95). 



of temperatures, where the emission in the lower excita- 
tion lines would have a relatively greater contribution from 
lower temperature gas. We repeat the caution that these 
ratios involve fluxes in the S(3) line which in rows 96, 97 
and 98 is blended with [Si VI], and it is these points which 
lie farthest from the LTE curve. The ratios of some para- 
hydrogen lines, viz. v=l-0 S(4)/S(2) and S(2)/S(0), also de- 
viate from the LTE prediction, thereby attesting to the ef- 
fect being real. These deviations indicate that the lines are 
either produced in low density gas (below the relevant crit- 
ical densities of ~ 10 5 cm -3 ) or that they are excited by a 
non-thermal mechanism, or perhaps both. Non-thermal ex- 
citation by suprathermal electrons liberated by hard X-ray 
photoionization within the cold, highly molecular core of the 
cloud was implicated by Jaffe & Bremer (1997) on the basis 
of the high v=l-0 S(l)/Ha ratio, and we consider it in detail 
below. 



4.2 Comparison with an X-ray dissociation region 
model 

Cold clouds deposited by cooling flows and immersed in the 
radiation field of the surrounding ICM, can produce molecu- 
lar line emission in two distinct regions. Part of the emission 
can be generated via thermal excitation immediately beyond 
the ionized skin of the cloud, peaking where the molecular 
fraction first reaches ~ 10~ 3 and the temperature ~ 2000 K. 
At greater depths, the temperature drops rapidly towards 
that of the microwave background and the gas becomes 
highly molecular (Ferland, Fabian & Johnstone 1994). At 
these depths, primary and secondary electrons generated in 
situ by hard X-ray photons can excite appreciable H2 emis- 
sion (this contrasts with the situation in a highly-ionized gas 
where primary photoelectrons rapidly thermalize through 
collisions with other electrons). Only by including this lat- 
ter contribution to the v=l-0 S(l) emission were Jaffe & 
Bremer able to reconcile the high observed v=l-0 S(1)/Hq 
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ratio with the calculations of Ferland, Fabian & Johnstone. 
The most detailed calculations of the emission from X-ray- 
irradiated molecular gas are those by Maloney, Hollenbach 
& Tielens (1996) (hereafter MHT) and we proceed to apply 
them to Cygnus A. 

As discussed by MHT, the dominant parameter control- 
ling the physical conditions is Hx/n, the ratio of the X-ray 
energy deposition rate per particle, Hx, to the gas density, 
n. The former is the integral over energy of F(E)a pa ., where 
F(E) is the local photon energy flux per unit energy interval, 
and <7 P a the photoelectric cross-section per hydrogen atom. 
For an X-ray source with a power-law spectrum of photon 
index F = 2 and 1-100 keV luminosity Lx = 10 44 L44 erg s _1 , 
located 100i'2 pc from the gas cloud, MHT give the following 
expression for Hx' 

Hx ~ 7 x HT 22 L 4 4 r 2 ~ 2 N^ 1 ergs" 1 (1) 

where iV att = 10 22 A r 22 cm 2 is the hydrogen column density 
attentuating the X-ray flux; they assume a column of at 
least JVatt = 10 21 cm -2 in order to exclude the normal UV 
photon-dominated region at the cloud surface. Figs. 6a and 
6b of MHT show the surface brightness, Q(H x /n), (from 
a cloud of column density 10 22 cm -2 ) in the [Fell] Al. 644, 
Br7, v=l-0 S(l) and v=2-l S(l) lines as functions of Hx/n 
at the cloud centre, for densities n = 10 3 and 10 5 cm -3 . 
These figures can be used to compute the total luminosity, 
Li, in line i from an area A c \ of cloud- face exposed to the 
continuum source, emitted in the column density interval 
(Ni, 22 -iV2,22)10 22 cm" 2 : 

Li — A c i / " $ti(H x /n) dN 2 2 (2) 

J jV 122 

using eqn. (1) to evaluate Hx/n as a function of N 2 2- Thus 
eqn. (2) can be used in conjunction with the observed v=l- 
S(l) luminosity in each row to derive A c i as a function of 
N-2,22 and cloud density, and thence the total mass of gas in 
the cloud, not merely that in warm H2. 

The radio sources in the CCGs considered by Jaffe & 
Bremer (1997) were weak enough to be negligible in compar- 
ison with the CF itself as an X-ray excitation source. This is 
not the case for Cygnus A, which houses a powerful obscured 
quasar: the Ginga spectrum of Ueno et al. (1994) exhibits an 
absorbed power-law component, with a photon index r ~ 2, 
a line of sight absorption Nu — 3.7 x 10 23 cm -2 and an 
unabsorbed 2-10 keV luminosity of 1.1 x 10 45 ergs _1 . The 
only circumstances under which this component would not 
be the dominant source of X-rays at the radii of inter- 
est would be if the line of sight from the quasar nucleus 
to the extended gas were obscured by the above column 
density or more, whilst the X-rays from the cooling flow 
were unattenuated. The off-nuclear increases in the v=l- 
S(l)/[FeII]A1.644 and v=l-0 S(l)/Bry ratios (Fig. g) are, 
with reference to Figs. 6a and 6b of MHT, consistent with 
the expected off- nuclear decrease in Hx/n from a nuclear ex- 
citation source. In fact, since only the narrow component of 
[FeII]A1.644 is produced in the same region as v=l-0 S(l), 
the v=l-0 S(l)/[FeII]A1.644 ratio relevant in this context 
increases more sharply than Fig. li indicates. 



4-2.1 The mass of off-nucleus v=l-0 S(l)-emitting clouds 

Following the procedure outlined above, we compute the 
total mass of gas in the v=l-0 S(l)-emitting clouds (ion- 
ized+neutral+molecular) in each off-nucleus pixel. We as- 
sume that all of the emission is from gas at the projected 
radius of the pixel, which strictly speaking means that the 
calculated masses are lower limits on the true values pro- 
jected on to each pixel. We perform the calculation for two 
limiting cases: firstly, where the gas is exposed to the unob- 
scured quasar X-ray spectrum, obtained by extrapolation to 
~ 100 keV of the power-law observed by Ueno et al. (1994); 
secondly, for the case where the X-ray spectrum is absorbed 
by 3.7 x 10 23 cm -2 on nucleus, in a circumnuclear torus. 
For each of the two densities modelled by MHT, we in- 
tegrate eqn. (2) from N\ t m = 0.1 to A^,22 =1, 10 and 
100; in some cases the minimum value of Hx/n plotted by 
MHT (10" 29 ' 4 ergs _1 cm 3 ) is reached before ^2,22 = 100, 
in which case a point is not shown for this column. The 
results for the case where the gas sees an unobscured spec- 
trum, shown Fig. ^a, demonstrate that the inferred mass is 
relatively insensitive to the assumed gas density. There is a 
stronger dependence on the stopping column density, espe- 
cially close to the nucleus, where the bulk of the emission is 
produced near the thermal-excitation peak within JV2,22 = 1. 
Reducing the stopping column below N 2 ,22 = 1 (but still 
N 2i22 > 0.1), however, does not substantially reduce the im- 
plied gas masses, so the lowest points in each pixel of Fig. [?] 
may be considered approximate lower limits. Noting that 
each pixel of the detector covers a projected area of 1.7 kpc 2 
in Cygnus A, we use Fig. to place lower limits on the pro- 
jected gas mass as a function of radius: in rows 95 and 97 
(i.e. within 1 kpc of the nucleus) E > 6 x 10 s M Q kpc -2 (this 
does not include the spatially- unresolved nuclear emission) , 
and in rows 91-94, 98 and 99 (i.e. within roughly 5 kpc of the 
nucleus) E > 3 x 10 6 M kpc~ 2 . If the distribution of this 
material is azimuthally symmetric, this implies that there 
are at least 2.3 x 10 s M of H 2 v=l-0 S(l)-emitting gas 
clouds within 5 kpc of the nucleus. We note that the higher 
mass points in Fig. ^ imply that much more gas could be 
accommodated without overproducing the flux in this line; 
and if the points in row 89 are azimuthally-representative, 
there could be a comparable mass at r > 5 kpc. Correction 
of the H2 v=l-0 S(l) fluxes for intrinsic extinction would 
increase the implied masses yet further. 

The fact that the surface mass density is not centrally 
concentrated in Fig. [7|a indicates that this gas bears no re- 
lation to the quasar nucleus, apart from the fact that its 
emission is excited by it. Rather, we favour an interpreta- 
tion in which this cold material is part of that deposited 
by the 25OM0yr _1 CF. In their deprojection analysis of 
the ROSAT PSPC data, Reynolds & Fabian (1996) found 
a mass deposition rate of ~ 35M0yr _1 within the small- 
est spatial bin (radius 20 kpc); on the assumption that the 
M(< r) cx r observed at larger radii continues to hold, 
M ~ 10 Mq yr _1 within 5 kpc. Thus, on a timescale of 10 9 
years, the CF would have deposited ~ 10 10 Mq within this 
radius. From Fig. we see that we can account for this 
quantity of material without violating constraints on the 
v=l-0 S(l) line emission only by integrating to a column 
density of 10 24 cm -2 . 

The results for the case where the X-ray spectrum is 
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Figure 7. The total mass of gas in the V=l-0 S(l)-emitting 
clouds, as functions of off-nucleus row number (the nucleus is 
in row 96), for densities of n = 10 3 cm -3 (open symbols) and 
10 5 cm -3 (filled symbols), under the following two scenarios: 
(a) when the emitting clouds are exposed to an unattenuated 
quasar spectrum, with the square, circular and triangular 
symbols denoting integration of eqn. (2) to stopping column 
densities of 7V H = 10 22 , 10 23 and 10 24 cm" 2 , respectively; (b) 
when the emitting clouds see a quasar spectrum absorbed 
by Nn = 3.7 x 10 23 cm -2 due to the nuclear torus; see text 
for details. Each pixel is 0.6 arcsec across. 



absorbed on nucleus are shown in Fig. [?]b. In this case, the 
clouds are exposed to a much weaker X-ray flux, and the 
emission is produced entirely in the non-thermal regime of 
Figs. 6a and 6b of MHT. We consider the case where the 
clouds are much less than 3.7 x 10 23 cm" 2 thick, which re- 
sults in an X-ray energy deposition rate which is relatively 
independent of depth within the cloud. The implied gas 
masses are thus independent of the actual cloud thickness. 
For the n = 10 5 cm -3 case, we calculate a total mass of 
1.6 x 10 10 Mq within a projected radius of 5kpc, which ac- 
counts for more than 10 9 yr of mass deposition from the 
10 Mq yr -1 CF within this radius. For the low density case, 
n = 10 3 cm -3 , the much higher mass would clearly violate 
kinematic constraints. 



4-2.2 The nuclear emission 

In addition, the nuclear v=l-0 S(l) flux can be used to 
place constraints on the putative circumnuclear torus, which 
may be responsible for the X-ray absorption towards the nu- 
cleus. According to Maloney (1996), this structure is likely 
to be ~ 100 pc in size: the detection of broad HI absorption 
(FWHM~ 270 km s" 1 ) towards the nuclear radio source 
(Conway & Blanco 1995) argues against a ~ kpc-scale torus; 
limits on the free-free absorption at 1.4 GHz and the effects 
of the non-thermal continuum on the spin temperature of 
the 21cm line rule out a pc-scale torus. We noted in sec- 
tion 3.3, however, that Imanishi & Ueno (2000) argued in 
favour of a < 10 pc-scale torus. 

Using a radius of 100 pc and integrating eqn. (2) to 
the X-ray-inferred column density of A r 2,22 = 37, we find 
a v=l-0 S(l) surface brightness of 2.3 erg cm" 2 s" 1 for n = 
10 5 cm -3 (for n = 10 3 cm -3 , the value of Hx/n in the cloud 
exceeds the range of the MHT models). In conjunction with 
the nuclear v=l-0 S(l) luminosity of 2.5 x 10 40 ergs -1 (cor- 
rected for light scattered into neighbouring rows) , we deduce 
a total gas mass of 3.3 x 1O 6 M . As seen from the central 
source, such material would cover 0.9 per cent of the sky. If 
the density of the torus is closer to n — 10 3 cm -3 , the molec- 
ular fraction will be much lower, leading to a correspond- 
ingly larger gas mass and covering factor. As discussed by 
Maloney (1996), the failure of Barvainis & Antonucci (1994) 
to detect CO absorption is most simply explained if the torus 
is wholly atomic. 

Alternatively, if the inner radius of the torus is < 10 pc, 
as suggested recently by Imanishi et al. (2000), Hx/n would 
be much higher and outside the range plotted by MHT. For 
densities of n = 10 3 and 10 5 cm -3 the surface brightness of 
the S(l) line would thus be much lower, leading to a higher 
gas mass and covering factor. Other excitation mechanisms 
may, however, be important, as discussed in section 4.3. 

Black (1998) has suggested that if the torus were mostly 
molecular, the H3" ion may be sufficiently abundant to pro- 
duce observable ro-vibrational lines in the K-band. He as- 
sumed the torus to lie at 300 pc from the nucleus, and calcu- 
lated the line intensities for two limiting cases, correspond- 
ing to kinetic temperatures of 100 and 1000 K. None of his 
predicted lines (which would have intensities in the range 
0.2-1.0 times that of Br7 in Fig. ^]) are seen in our spectrum, 
which also suggests that the torus is not highly molecular. 



4.3 Alternative excitation mechanisms 

We now consider two other excitation mechanisms for the 
H2 emission, namely UV radiation from young stars or 
the active nucleus itself, and shocks. In low density gas 
(rt <C 10 5 cm" 3 ), rovibrational H2 emission can be produced 
by direct photoexcitation of excited vibrational states, char- 
acterised by strong emission in the v=2-l S(l) line relative 
to the v=l-0 lines of lower excitation; in gas above the criti- 
cal density (~ 10 5 cm" 3 ), a sufficiently strong UV radiation 
source can heat the gas to a kinetic temperature of ~ 1000 K, 
leading to collisional excitation of the H2 lines (see the de- 
tailed calculations of Sternberg & Dalgarno 1989, hereafter 
SB89). 

We first consider the energetic requirements for UV ex- 
citation. The average surface brightness of the v=l-0 S(l) 
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line over rows 89-99 is 2.6 x 10 -5 erg cm~ 2 s" 1 sr^ 1 in the 
Cygnus A frame (values for individual pixels are within a 
factor of 4 of this value). From the models of SB89, the 
emission in this line represents no more than 2 per cent 
of the total 1-4/wn H2 line flux, so the overall H2 line 
surface brightness is at least 1.3 x 10~ 3 erg cm~ 2 s _1 sr _1 ; 
from Figs. 9 and 10 of SB89, this requires a density n > 
10 4 ' 7 cm -3 , and a 5-13.6 eV UV background which is a fac- 
tor of x si 10 3 higher than the standard galactic value of 
2.5 x 10" 4 erg cm -2 s _1 sr _1 (Draine 1978). The implied 5- 
13.6 eV luminosity is 1.4 x 10 11 L Q . 

Star-formation is often observed in the central galaxies 
of CF clusters (Allen 1995; Crawford et al. 1999), so we con- 
sider it as a possible source of exciting UV radiation. From 
HST WFPC2 imaging, Jackson et al. (1998) identified sev- 
eral blue condensations within the the spiral structure which 
is present within a 2 arcsec radius of the nucleus of Cygnus 
A, and deduced from their colours that they are plausibly 
due to star-formation within the last Gyr. Their implied lu- 
minosity of 3 x 10 7 L-0 is, however, insufficient to power the 
observed H2 emission. Furthermore, our north-south slit lies 
along the nuclear dust-lane and does not encompass any of 
these condensations (the finite seeing will of course scatter 
some of their light into the slit). On energetic grounds we 
cannot completely rule out the possibility that the H2 emis- 
sion is powered by heavily obscured star-formation along the 
dust-lane, since the monochromatic IRAS 60/im luminosity 
of Cygnus A, although very weak for its radio lobe power 
(Barthel & Arnaud 1996), is 5.25 x 10 11 L . 

The active nucleus is also a potential source of exciting 
UV photons. Imaging polarimetry by Tadhunter, Scarrott & 
Rolph (1990) and Ogle et al. (1997) suggests that there is 
an extended region which scatters light from a hidden nu- 
cleus; its morphology resembles that of the ionization cone 
seen in the HST image of Jackson et al. (1996), with a half 
opening angle of 55 degrees, oriented approximately along 
the radio axis. Shaw & Tadhunter (1994) performed long-slit 
spectroscopy to determine the spectral energy distribution 
(SED) of the featureless blue continuum (FBC) in the nu- 
clear regions, after removing the old galaxy component. In 
the region covered by our slit, they found a monochromatic 
surface brightness of 1.4 x 10 -5 erg cm -2 s _1 sr _1 A -1 at a 
rest-frame wavelength of 3620A, after correction for intrinsic 
extinction of E(B-V)=0.9 mag; the SED of the FBC steepens 
below 4200A, in a manner which could be consistent with the 
scattering of quasar light by dust. If Fx continues to be at 
least as steep as A -3 below 3620A, then its 5-13.6 eV surface 
brightness will exceed 0.35 erg cm -2 s _1 sr -1 , and thus be 
sufficient to power the observed H2 emission. Furthermore, 
since the FBC is scattered nuclear light, any H2 within the 
ionization cone will be exposed to a stronger UV radiation 
field (which could, however, be sufficient to destroy the H2 
altogether). 

Thus, on energetic grounds, scattered nuclear UV light 
or obscured star formation could potentially power the ob- 
served H2 emission. Both possibilities are, however, largely 
ruled out by the observed line ratios. As mentioned earlier, 
the v=2-l S(l)/v=l-0 S(l) flux ratio is an important diag- 
nostic for UV excitation. The v=2-l S(l) line is detected 
only in rows 94 and 96, with v=2-l S(l)/v=l-0 S(l) flux ra- 
tios of 0.23±0.06 and 0.13±0.03, respectively, for which the 
equivalent collisional (LTE) excitation temperatures, T12, 



are 3130 and 2380 K. For UV pumping in low density gas , 
T12 ~ 10 4 K. With the higher density gas (n > 10 4 ' 7 cm" 3 ) 
and intense UV radiation field (x ^ 10 3 ) required on en- 
ergetic grounds, T12 ~ 1000 K and the v=l-0 transitions 
would be collisionally-excited at the kinetic temperature of 
the gas (also ~ 1000 K), contrary to the results displayed in 
Fig. ^. We thus conclude that UV radiation does not play 
an important role in the excitation of the H2 emission. 

Finally we consider shocks as an excitation mechanism, 
following the method of Draine & Woods (1990). The frac- 
tion of the kinetic energy dissipated in a shock which is 
radiated in the v=l-0 S(l) line is e(v a ), where v a is the 
shock speed. For magnetised molecular clouds, e ~ 0.02 for 
30 < v B < 50kms _1 (Draine, Roberge & Dalgarno 1983), 
decreasing at lower shock speeds because the gas is not hot 
enough to emit in this line, and at higher shock speeds be- 
cause the H2 is dissociated. For the v=l-0 S(l) slit lumi- 
nosity of 5.7 x 10 40 ergs -1 , we thus require a power of ~ 
3x 10 42 erg s _1 to be dissipated in shocks. The kinetic energy 
of ~ 10 9 Mq of gas (equal to that which the CF would have 
deposited within the slit over 1 Gyr), moving with random 
velocities of ~ 50 kms -1 , is ~ 3 x 10 55 erg; it could thus sup- 
port this rate of dissipation for only ~3x 10 yr. Since the 
mass of warm molecular hydrogen is only 2.4 x 10 5 Mq (see 
section 4.1), only ~ 10~ 4 of the cold material from the CF 
would be emitting at any one time. As stated in Jaffe & Bre- 
mer (1997), assuming that it takes lyr for shocked material 
to return to quiescence, the material must be shocked every 
10 3 — 10 4 yr. Based on the shortness of these timescales and 
on the mismatch between the required random cloud veloci- 
ties for efficient H2 emission (~ 50kms _1 ) and the observed 
line widths, we do not consider shocks to be important in 
this context. Hard X-ray heating remains the most viable 
excitation mechanism. 



5 SUMMARY AND CONCLUSIONS 

We have presented J, H and K-band UKIRT CGS4 spec- 
troscopy of Cygnus A, as part of a programme to study the 
H2 emission in a sample of cooling flow central cluster galax- 
ies. Spanning a rest-frame wavelength range of 1.0-2.4^tm, 
the spectra include the S(0)-S(5) lines of the H2 v=l-0 se- 
ries, Paa,7, Br7,<5 and [Fell] AA1. 258, 1.644. The lines are 
spatially extended by up to 6 kpc from nucleus. Kinematic 
differences between them indicate that the three types of 
line emission (H2, H recombination and [Fell]) are not pro- 
duced entirely in the same gas clouds. In particular, much 
of the [FeII]A1.644 line originates in a broad component 
(FWHM ~ 1040kms _1 ), perhaps due to interaction with 
the radio source, but an unresolved core also contributes 
on nucleus; the H recombination lines and those of H2 are 
succesively narrower. Foreground screen extinction models 
imply Av=9.6 (from Br7/Paa), 5.7 (from Paa/Pa7) and 12 
mag (from the [Fell] lines), substantially greater than the 
extinctions deduced by previous workers using optical lines; 
comparable extinctions are also deduced in some off-nucleus 
rows. 

The relative intensities of various H2 v=l-0 S lines devi- 
ate from the expectations for thermal (collisional) excitation 
in gas at a single kinetic temperature. Hard X-rays from the 
obscured quasar nucleus are found to be the most plausi- 
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ble excitation source for the H2 emission; we thus applied 
the results of MHT on molecular line emission from X-ray 
dissociation regions to deduce the projected mass density 
(E) in v=l-0 S(l)-emitting clouds, for gas densities of 10 3 
and 10 5 cm -3 and stopping column densities between 10 22 
and 10 24 cm -2 . If the observed variation of £ along the slit 
is azimuthally representative, it implies a mass of at least 
2.3 x 10 8 Mq within 5kpc of the nucleus. The near uniform 
variation of E with radius strongly suggests that this ma- 
terial represents part of that deposited by the cooling flow 
(for which M ~ 10M Q yr _1 inside this radius) rendered 
visible in this manner by the 'searchlight' of the powerful 
quasar. Moreover, if much of the X-ray absorption to the 
point source seen by Ueno et al. (1994) actually arises in a 
< lOOpc-scale circumnuclear torus, the implied mass in cold 
clouds is much larger, ~ 10 10 M , plausibly accounting for 
~ 10 9 yr of mass deposition from the cooling flow. 

In a subsequent paper we shall apply these techniques 
to the larger sample of CF CCGs mentioned in section 1. 
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